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Abstract 

We explore the interactions of two strangeness -3 baryons in multiple spin channels with lattice 
QCD. This system provides an ideal laboratory for exploring the interactions of multi-baryon 
systems with minimal dependence on light quark masses. Model calculations of the two-O - system 
in two previous works have obtained conflicting results, which can be resolved by lattice QCD. 
The lattice calculations are performed using two different volumes with L ~ 2.5 and 3.9 fm on 
anisotropic clover lattices at m n ~ 390 MeV with a lattice spacing of a s ~ 0.123 fm in the spatial 
direction and at ~ a s /3.5 in the temporal direction. Using multiple interpolating operators from 
a non-displaced source, we present scattering information for two ground state f2~ baryons in 
both the S=0 and S=2 channels. For S=0, fccot<5 is extracted at two volumes, which lead to an 
extrapolated scattering length of a^Oo = 0.16 ± 0.22 fm, indicating a weakly repulsive interaction. 
Additionally, for S=2, two separate highly repulsive states are observed. We also present results on 
the interactions of the excited strangeness —3, spin-^ states with the ground spin-| states for the 
spin-1 and spin-2 channels. Results for these interactions are consistent with attractive behavior. 

PACS numbers: 12.38.Gc, 12.39.Fe 
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I. INTRODUCTION 



Lattice QCD calculations have advanced to the point that scattering phenomena for 
multi-hadron systems can be reliably calculated from first principles. These calculations, 
performed through the analysis of two or more hadrons in a finite volume, allow for phase 
shifts and potential bound states to be studied non-pertubativelyp], 0]. The majority of the 
focus of these calculations has been to explore baryon-baryon and meson-meson systems, 
where in the latter the scattering length for the 1=2 irir scattering system has been calculated 
to within a few percent [SMll . Additionally, high precision calculations of bound states for 
the two lambda system|10l. 111] , as well as the deuteron and SS system 12J were performed 
recently. In this work, we explore a different hyperon-hyperon system, namely the two 
strangeness -3 baryon (QQ) system, where we present the scattering results for both ground 
and excited states. 

While lattice QCD calculations of excited states using cubic irrep sources are well es- 
tablished, these lattice techniques have only been applied to mesonic scattering^ and have 
yet to be applied to two baryon scattering. Single hadron excited states have seen a great 



deal of attention from the lattice community for both mesons [13l. Il4j and baryons |l5Hl 7 



where many states are consistent with their experimental counterparts. Recent advances 
include calculations of the isoscalar meson spectrum |l4fl with the use of the latest algo- 
rithmic methods for calculating disconnected diagrams [18| and the employment of GPU 



technology [19[. We extend this approach of utilizing multiple embeddings of lattice ir- 



reducible representations |20| to the two baryon system, with the ultimate goal being to 
extract higher partial wave interactions of nucleons from fundamental lattice calculations. 
As a starting point, we apply these techniques to the QQ system in a relative s-wave state. 

The lattice study of the flfl system is of interest for several reasons. Like most hyperon- 
hyperon systems, the interactions between two or more Q~ baryons are poorly understood 
experimentally due to their large mass and relatively short lifetime. To this end, lattice 
QCD calculations can predict phenomena in these systems and pinpoint signals for heavy 
ion scattering experiments, such as STAR or ALICE. The QQ system has not received as 
much theoretical attention as its lighter hyperon counterparts, such as the H-dibaryonj2l| 
and single A hypernuclei j2~2h 24l| . However, within the last decade, this system was studied in 
the context of the chiral quark model[25[, where it was found to prefer a bound ground state 
with a binding energy on the order of 100 MeV. A conflicting analysis(26| using the quark 
dislocation model finds the system to be weakly repulsive. Additionally, the interactions of 
the QQ system, along with the interactions between the Q~ and other baryons, may prove 
to be relevant in dense systems several times nuclear density, such as the core of a neutron 
star 27J ■ These interactions may also play a role in the core of a supernova, ultimately 
determining whether the system becomes a neutron star or a black hole. 

Another attractive aspect of studying the QQ system on the lattice is the fact that the 
system is believed to primarily depend on the physical strange quark mass as opposed to the 
unphysically large light quark masses with which these lattice calculations are performed. 
This assertion is found to be true for the f2~ and several of its excited states in Ref. Il6| 



and, consequently, the f2 mass is often used to set the lattice scale 28] . Thus, unlike most 



nuclear calculations involving light quarks, calculation of two-f2~ systems at the physical 
point should rely less on chiral extrapolations. Additionally, the inversions involving only 
the strange quark are less computationally expensive and the resulting signal involving 
strange baryons is cleaner. For these reasons, the multi-f2~ system is the ideal laboratory 
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for understanding nuclear interactions on the lattice directly as it not only could provide 



insights into two or more nucleon interactions [12|, I29H31 



opportunity to study a host of nuclear interactions, 



but it also provides a unique 
such as the tensor forces and s-wave 



three-baryon forces |32l. 133 



This work is organized as follows: In Sect. [TT1 the basic properties of the QQ system are 
mapped out along with differences from the two nucleon systems. In Sect. IIIH the QQ system 
in a finite box is explored, conventions are defined, and the methods of multiple embeddings 
are discussed. In Sect. IIVt the calculation details and analysis methods are explained, with 
the lattice results presented in Sect. [V] and Sect. [VH Finally, in Sect. IVIII the scattering 
results are derived, with a conclusion in Sect. IVIII1 



II. TWO BARYON SYSTEM IN INFINITE VOLUME 

As in the case of two nucleons in the isospin limit, the channels with which two-f2~ 
baryons can interact are restricted by the Pauli exclusion principle. Each Q~ contains three 
valence strange quarks and is spin-3/2 in its ground state. Pauli statistics dictate that 
the two Q~ wavefunction must be antisymmetric. Where the QQ system differs from two 
nucleons is each Q~ is spin-3/2 as opposed to spin- 1/2 and there is no isospin wavefunction. 
In general, two spin-3/2 particles can exist in a total spin S — 0, 1, 2, or 3 state. However, 
the additional condition for an anti-symmetric wavefunction leads to the result that two Q~ 
baryons in an s-wave state can only have two non-trivial spin channels: S=0 and S=2. 

As mentioned previously, there have been two model calculations of the QQ system in the 



S=0 channel. The first model calculation in Ref. [26[ explores potential di-baryon systems 
via the quark dislocation and color screening model. In this model, the quark model is 
generalized in several notable ways. Namely, the color screening/string tension interaction 
is included in the Q-Q potential along with a delocalized quark orbit in the wavefunction. 
The authors find good agreement with the experimental N-N system. For the QQ system, 
the authors find a weakly repulsive interaction (26| 

AE nn = E nn - 2M n = 43 ± 18 MeV. (Quark Disloc. /Color-screen Model) (1) 

Here, the positive AE value implies the theory is not bound and likely weakly repulsive 
(near threshold). 

The calculation in Ref. [25| explores the QQ system in the chiral SU(3) quark model. 
In essence, the chiral SU(3) quark model generalizes the quark model, consisting of one- 
gluon exchange and a confining potential, to an SU (3) sigma model in order to account for 
non-perturbative effects that affect the constituent quark mass. The resulting Hamiltonian 
from the confining potential in this set up (whose 17 free parameters are determined from 
experiment) allows for two baryons systems to be studied by solving the resonating group 
method equations. The authors fi nd g ood agreement with experimental N-N and Y-N phase 
shifts, and find for the system [251] 

AE nn = E nQ - 2M n = -116 MeV. (SU(3) Chiral Quark Model) (2) 

The depth of this bound state is significant and could be detected in heavy ion experiments, 



as detailed in Ref. [25]. Also, a binding energy of this magnitude would easily be resolved 



in state-of-the-art lattice QCD calculations by multiple standard deviations. Reference 



251 ] also points out potential issues with the quark derealization model due to non uniform 
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FIG. 1: Leading diagram involving pion exchange for the 00 system. 



confinement potentials [34J. Ultimately, one would prefer a first-principle, model-independent 
lattice QCD calculation to address this debate. 

An over-arching goal of lattice studies of multi-baryon systems is to explore the connec- 
tions between lattice calculations and parameters of many-body effective field theory. This 
has been a primary goal for light baryons, but due to unphysically large quark masses this 
connection has proved difficult to achieve. In the case of the multi-f2~ system, the light quark 
mass dependence is expected to play a minimal role, as the leading order interaction involv- 



ing pions is given by processes involving pair-produced two pion vertices [35[], as depicted in 



Fig. [TJ Consequently, lattice calculations with unphysical light-quark masses should provide 
'near-physical' results in the f2-only sector. To date, there has not been much development 
in multi-f2~ EFT due to the fact that low-energy physics of weakly decaying Q~ baryons 
is difficult to probe physically. However, with lattice QCD calculations, acquiring parame- 
ters for a meson-less EFT of multiple Q~ baryons should be possible. Recent workj36| has 
touched on this subject within the context of dark matter, and in a forthcoming paper we 
develop in detail the two flavor EFT for the two-f2 system. 

In this work, we examine both ground state interactions and interactions of excited states 
in the s-wave two Q~ system at a single pion mass of 390 MeV with two separate volumes. 
Ultimately, we intend to probe tensor interactions by projecting the initial and final state 
baryons to higher partial waves. Future calculations will quantify to what extent the claim 
of minimal dependence on the light quark mass is valid. Once understood, a new gateway 
to understanding interactions and systematics between EFT and lattice QCD calculations 
can be probed in a way currently inaccessible to light baryon systems. 



III. TWO SYSTEM IN A FINITE BOX 

There are several well-known complications when studying scattering calculations on the 
lattice. First, discretized lattice calculations with periodic boundary conditions no longer 
preserve the 0(3) rotational symmetry observed in the continuum, but rather preserves an 
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octahedral subgroup. Thus, in order to extract information about the continuum S=0 and 
S=2 states, the appropriate irreducible representations of the octahedral group must first be 
understood. Second, lattice calculations are performed in Euclidean space, where the usual 
LSZ formalism only holds at kinematic threshold [37]. As a result, phase shifts have to be 
extracted by analyzing energy shifts of hadrons in finite volume. For two hadrons (A and 
B) in a finite volume, the energies associated with the four-point correlation function are 
given by 



E = Jk 2 + m 2 A + Jk 2 + m% = AE + m A + m B , (3) 

where AE is the energy of interaction and k is the associated momenta in the center of mass 
frame. As was shown in Refs. [l|, 0], this interaction momenta at a given volume can be 
related to scattering phase shifts by 



kcot5(k) = —S[ — , (4) 




where the function S is given by the regularized three-dimensional zeta-function 

|j|<A 

s w = £ - 47rA - ( 5 ) 

The value of k cot 5(k) extracted then has the normal effective range expansion given by 

k cot 5(k) = -- + ^rk 2 + ... (6) 

and given multiple values of k cot 5(k) at multiple values of k 2 (obtained through lattice 
calculations of the same system at differing lattice volumes or through the use of boosted 
systems) one can make an extraction of the specific scattering parameters a, r, and so on. It 
is important to note that the relation in Eq. H] holds for elastic scattering processes and no 
longer holds when the energy of interaction exceeds twice the pion mass, while the relation 
in Eq. [6] is only valid below the t-channel cut. These facts often limit the number of excited 
scattering states that can be extracted from the lattice at a fixed volume. The best way 
to understand this is through the realization that two non-interacting hadrons in a finite 
volume have quantized non-relativist ic energy levels given by 4irn/(mL 2 ) for select integer 
values of n. At a fixed volume, this reduces the access to excited scattering states for light 
hadrons, but allows for more excited scattering states to be explored for two Q~ baryons. 
Thus, for extracting excited behavior between multiple hadrons, including higher partial 
wave and tensor interactions, the two Q~ system is superior. However, as will be discussed 
in more detail in the following sections, having more excited states accessible leads to more 
excited state contamination when trying to extract a given state. 

In order to extract information about the Q~ , one must first calculate using an operator 
that has some (preferably large) overlap with the ground state. Systems with definite integer 
spin modes in the continuum limit have dominant overlaps with different lattice irreducible 
representations (irreps), labelled by Ai, A 2 , E, T 1; and T 2 . The same can be said for fermionic 
modes which can map on to the lattice irreps Gi, G2, or H. Further, each representation T 
will have an associated parity, which we label as for postive or negative parity. Table |T] 
enumerates the different lattice irreps used in this work and provides their dominant spin 
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TABLE I: Lattice irreps T used in this work and dominant overlap with spin J in the infinite 
volume limit (L — > oo). 



r (L ^ oo) 


J (L = oo) 


At 





T{ 


1 


E+ 


2 


r+ 


2 


Gf 
H+ 


l 

2 
3 
2 



content in the infinite volume limit. As seen from this table, the spin-3/2 Q~ particle is 
represented by the H + irrep in finite cubic volumes. 

The interpolating operator representing the Q~ baryon is given byj20| 

= £abc S a S /3 S y> CO 

where a, b, and c are color indices and a, 0, and 7 are spinor indices. Appropriate linear com- 
binations of the spinor indices will produce Q~ interpola ting operators with definite lattice 
symmetry. For non-displaced sources and following Ref. [20], there are two representations, 
or embeddings, of the Q~ particle in the H + irrep given by 



H + : 



Interpolating operator 


\J,Jz) 




Interpolating operator 


\J,Jz) 




|3/2,3/2) 




\/3 ^133 


13/2,3/2) 




|3/2,l/2) 


2 H + : 


2O134 + f^233 


13/2,1/2) 


^122 


|3/2,-l/2) 




^144 + 21^234 


|3/2,-l/2) 


^222 


|3/2,-3/2) 




\/3 f^244 


|3/2,-3/2) 



For a given embedding ( 1 H + or 2 H + ), each infinite volume \J, J 2 )-state and its corresponding 
source are given. It is important to note that each state within a given embedding is 
orthogonal to the other states in that embedding after averaging over configurations. Thus, 
each embedding can lead to four statistically independent calculations of the Q~ two point 
function. In the non-relativistic limit, the first embedding maps onto the upper two spinor 
components in the Dirac-Pauli basis, while the second maps onto the lower two components. 
As such, one expects larger overlap with ground state systems when dealing with the first 
embedding. It is also important to note that contracting the same state between two different 
embeddings is statistically dependent. 

The s-wave states of the two Q~ system can be formed from a tensor product of the 
ground state lattice irreps. The lattice irrep that has overlap with the S=0 state is the 
Af irrep. In terms of the spin-3/2 states of the fl~ ground state, there is only one linear 
combination that leads to Af and it is given byj20| 

At(S = 0) : - ( HsH_s - H_sHs + H_i Hi - H L H_i) (8) 

- 1 - 2 V 2 2 22 22 22/ 

where the subscript indicates the z-component of the spin and the if's can be in either 
embedding. 
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For S=2, there are two lattice irreps that have overlap with the ground state, namely E + 
and T 2 + , and as a result both lattice irreps can be used as independent calculations of the 
S=2 QQ scattering state. Additionally, unlike the A± lattice irrep, both E + and T 2 + can be 
formed by multiple linear combinations of the two H + states, each of which are statistically 
independent determinations of the lattice irrep. For E + , there are two linear combinations 
given by 

E + (S = 2) : 
1 : 



H_ 3 Hs 

2 2 



2 : 



,HiH a 

2 V 2 ~2 

— I Hi Hi — Hi Ha 

2 V 2 2 2 2 



H_iHi 4- HiH_i 

2 2 2 2 



H _ 3 H_ i 

2 2 



H_iH_ 3 

2 2 



(9) 



and for T 2 + , there are three linear combinations given by 

T+(S = 2) : 
1 : 



-^=(h 3 h_i 

V2\ 5 2 



H_iHz 

2 2 



2 
3 



-[Hi Hi 

2 V 2 2 
-^=(HlH_i 

y/2\ 5 2 



#1#3 

2 2 



H_iH_i+ H_iH_i 

2 2 2 2 



H^sHi 

2 2 



(10) 



It is worth noting that all linear combinations for E + and T 2 + above will yield zero if each 
H + irrep is in the first embedding (the l.h.s. of Eq. flH])). Thus, to extract S=2 state, one 
minimally needs one H + irrep in the first embedding and the other in the second embedding, 
which naturally leads to a higher level of excited state contamination. 

A good check worth pointing out is that the S=l system should be trivial. More specif- 
ically, the linear combinations that form the Ti irrep should be zero due to parity re- 



strictions and anti-symmetry. Following Ref. [20j, one such linear combination is given 
by 3(H 3 / 2 H_i/ 2 + H-1/2H3/2) —4(^1/2-^1/2)- Since each H represents a source of three iden- 
tical strange quarks, the exchange of any two H terms will lead to an overall minus sign. 
Thus, the first two terms will exactly cancel with each other and the last term can only be 
zero. Thus, this linear combination will yield zero for the two Q~~ system. 

In addition to the H + irrep that represents the Vt~ ground state, one can additionally 
explore excited states that correspond to spin-1/2 in the continuum limit. The lattice irrep 
associated with this excited mode is Gf, which is given in terms of the Q~ operators 



Gt 



Interpolating operator 




J,Jz) 


^134 — ^233 
^144 — f^234 




11/2,1/2) 
1/2,-1/2) 



Since the H + baryon and the G^ baryon are not identical particles, scattering between these 
states can take all the spin values allowed by the addition of angular momenta (S=l,2). The 
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TABLE II: Gauge configuration details 



Size 


mi 


m s 


m n [MeV] 


m n L 


Configs x Meas/Config 


20 3 x 256 


-0.840 


-0.743 


~ 390 MeV 


4.855 


1155 x 50 


32 3 x 256 


-0.840 


-0.743 


~ 390 MeV 


7.74 


465 x 25 



s-wave S=l contribution (in irrep T x + ) is given by 

T+(S = 1) : 



1 : — ( GiHi — 3G_iUa 

2 \ 2 2 2 2 

2 : — —( G iUi — G_iHi 

V2\ 22 22 

3: -(3GiH_3-G_iH_i) (11) 



and the s-wave S=2 contribution (in irrep E + ) is given by 

E + {S = 2) : 



— ;= ( G 1H1 + G _i Hi 

2 2 2 2 



v/2 
1 

7! 



(12) 



Finally, it is possible to couple two G\ excited fl particles to form states of definite Ai 
symmetry, corresponding to an s-wave, S=0 system in the infinite volume limit, using the 
coupling coefficients given in Ref . [20| . However, as we point out below, because of the high- 
energy levels associated with this system, coupled with limited statistics, we were not able 
to extract any statistically meaningful information from this system on the current lattices 
that were used in this work. 



IV. LATTICE DETAILS 
A. Configurations 

Our calculations were performed on anisotropic Wilson lattices generated on the uBGL 



machine at LLNL using the tuning parameters defined in Ref. [28J. The primary ensemble 
used was a 20 3 x 256 with m n w 390 MeV, a s ~ .1227 fm, and a s /a t ~ 3.5 (see Ref. [12| for 
more details of the anisotropic parameters). The spatial extent of these lattices is L ~ 2.5 
fm. On the 20 3 x 256 ensemble, 50 propagators with random sources were calculated on 
every 5 trajectories, where every measurement required one propagator. The propagators 
were generated on the Edge GPU cluster at LLNL. For the S=0 two-Jl system, we also 
performed measurements on m v ~ 390 MeV lattices at a larger volume (32 3 x 256) m an 
attempt to quantify volume effects. Here 25 measurements per configuration, blocked every 
4 trajectories, were made. Table [TT] gives details of the configurations and measurements 
used for this work. 
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B. Contractions 



Due to the three degenerate valance quarks, the f2~ interpolating operator has several 
symmetries worth noting. The first symmetry is that the spinor indices can be interchanged 
freely, 

f2 Q/ 3 7 = £lp ai = VLp ia = ■ ■ ■ . (13) 

This fact is the result of exchanging two quarks leads to both a minus sign from permuting 
Grassman number and a minus sign from exchanging indices in the epsilon tensor, which 
cancel. As a result, all Wick contractions will have the same relative sign. For the two-Q 
system, it can be shown using these symmetries that all contractions fall into two distinct 
forms: 'direct' and 'exchange', as shown in Fig. |2j Taking advantage of these symmetries 



{a 1 a 'l C a » a' | 

P ' P'> <W %y < P ' P'> "„ 




r - 1 r ^ i 

I CO > CO J I CO * CO j 



(a) (b) 

FIG. 2: Quark contractions types used in the calculation including the (a) direct contributions and 
the (b) exchange contributions. 

drastically reduces the 6!=720 possible contractions of two-f2 system by an order of mag- 
nitude and therefore reduces both computational cost due to matrix multiplication and 
required memory to hold the matrices. 

Before performing the contractions, we first calculate the propagator using a Gaussian 
smeared source. Upon inverting on that source, we Gaussian smear the sink (referred to 
hereafter as shell-shell, or SS measurements) or perform no smearing to the sink (referred 
to hereafter as shell-point, or SP measurements). While the SS measurements are expected 
to give the best overlap with the ground state spectrum, by utilizing the SS and SP mea- 
surements in combination, one can largely eliminate the contribution from the first excited 
state and extract a more robust ground state signal that dominates at earlier Euclidean 



time|38l. |39[. One can further enhance the ground state signal by making use of different 
combinations of the embeddings 1 discussed in Sect. Ill II Using the propagator, we form the 
relevant irrep "blocks" for a given embedding, and we further suppress the excited states 
by projecting the momentum of both individual irrep blocks to zero independently, which 
will result in a correlation function where excited states with nonzero back-to-back momen- 



tum have been removed [301] . Finally, these irrep blocks are contracted and the correlation 



functions of interest are obtained. 



Previous studies 



15-17 



28 



40( have used a matrix of correlation functions from different embeddings to 
extract the low-lying spectrum of several baryons, up to the first few excited states. As we are primarily 
concerned with the ground state baryons and their interactions these techniques are not used here. 
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C. Analysis Details 



The measurements are blocked by configuration and the resulting ensemble is boot- 
strapped, with each bootstrap measurement and the final bootstrap ensemble being the 
same size as the original ensemble. For each embedding correlation function, the SS and SP 
measurements are put through a matrix-Prony algorithm as detailed in Ref. 38j. Specifi- 
cally (following the notation and derivation of Ref. 38] ), the correlation function recursion 
relation 



My r {t + tj)-Vy T (t) = 



has as a solution 



M 



'r+t w 



yr(t + tj)yr(t) rj 



t=T 



V 



T+tw 



J2 Vr(t)Vr{t)' 



t=T 



(14) 



(15) 



for the vector of correlation functions yr(t), of irrep and embedding type T. Here, the 
window of timeslices from r to r + t\y is the set of values of the correlation function over 
which the outer product is taken. This window must include enough information to make 
the resulting matrix full-rank and invertible, with subsequent timeslices helping to reduce 
statistical noise to some extent j38|. The choice of tj will increase the "lever-arm" that 
the matrix-Prony rotation provides, further mitigating statistical noise but at the price of 
increasing systematic fluctuations in the correlation function. As such this integer quantity 
should typically be chosen to be small. The eigenvectors qr of the matrix V~ l M, are defined 
by 



V~ l Mq x 



A 



(16) 



where the eigenvalues Ar,j are placed in ascending order. This will return a correlation 
function gr,o that has an enhanced ground state contribution. Using the rotation matrix 
defined by qr on each of the bootstrap measurements for type T will lead to an effective 
mass plot of the function 



M eff) r(t) 



1 



9r,o(f) 



(17) 



with a longer and more robust plateau region. For each of the effective mass plots below, the 
values tw = 10 and tj = 2 have been used, with different rotation points r for each type T. 
A fully correlated x 2 minimizing fit is then performed in the plateau region on the ensemble 
of bootstrapped effective mass data to extract the ground state energy with a statistical 
error. A systematic error from the fit window choice is obtained by modifying the endpoints 
of the fit window ±2 timeslices and taking one-half the maximum minus the minimum of 
those fit values. The fit values are displayed along with the \ 2 /dof for the fit and the Q (or 
quality of fit) value, which is the integrated probability distribution of x 2 from the observed 
fit x 2 to infinity. 

With this (and other black-box) methods there is serious concern that a false plateau 
may be recovered, as competing overlap factors under the rotation may produce a temporary 
cancelation in the plateau region that mimics the behavior of a real plateau. This behavior 
may be particularly problematic for methods with a large number of correlation functions in 



10 



the vector yr, as the opportunities for cancelations of the wrong form increase. To remove 
this complication, in this work all plateaus are initially identified in the effective mass fit 
of the SS data for each irrep and embedding T. This correlation function is manifestly 
positive (up to an overall sign) and so does not suffer from possible overlap cancelations. 
The matrix-Prony rotation point r is then chosen such that gr,o returns an effective mass 
plateau fit value within 0.5a of the SS plateau fit value. The matrix-Prony result will then 
return a value that is statistically the same as the real plateau, but with a significantly 
improved signal-to-noise ratio. 



V. STRANGENESS -3 SYSTEMS 



A. 



§: H + Irrep 



As discussed in Sect. IIIH the strangeness -3 H + irrep has a dominant overlap with a 
spin-3/2 particle, whose ground state is the Q~ particle. This irrep has two embeddings, 
making for a total of four possible source/sink embedding combinations. We will use the 
notation for correlation functions with source embedding type i and sink embedding 
type j. In Fig. El these four combinations are plotted for both the pure SS data and for data 
that has undergone a matrix-Prony rotation at timeslice r = 15 for the 20 3 x 256 lattices. 
While the largest overlap is observed in the combination, it is apparent from the figure 
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FIG. 3: (Color online) Effective mass plots for the four different H + embedding combinations 
calculated using the 20 3 x 256 lattices, using (a) the pure SS data and (b) a matrix-Prony rotation 
about timeslice r = 15. The four embedding combinations shown are (blue circles), H^ 2 (red 
squares), (brown diamonds), and (green triangles). 



that the highest embedding combination, H^-, has a significant amount of excited state 
contamination. Furthermore, this embedding combination does not have enough overlap 
with the ground state for a signal to appear before the onset of the baryonic noise around 
timeslice t = 40. This lack of overlap with the ground state persists even following the 
matrix-Prony rotation, indicating that this higher embedding combination has an overlap 
with the ground state that cannot be resolved with the statistics available for this calculation. 
A similar situation is observed for the 32 3 x 256 lattices. 
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Given the high number of measurements made for this calculation, this finding throws 
into doubt the utility of H22 for any method attempting a better extraction of the Q~ ground 
state. Specifically worrisome is the possibility that the contribution of even small amounts 
of this correlator to an effective mass plateau may lead to an inaccurately high ground state 
energy. As such, this embedding combination is removed from consideration in the following 
H + discussion. The remaining three embedding combinations have significant overlap with 
the ground state, and the effective mass plot for the sum of the bootstrap ensembles of the 
three lowest embedding combinations is shown in Fig. H] for both the 20 3 x 256 and the 
32 3 x 256 lattices. The results of the fits to each data set are shown in Table IIHI 




FIG. 4: (Color online) Effective mass plots for the H + (S=|) Vt~ baryon calculated using (a) 
20 3 x 256 and (b) 32 3 x 256 lattices. The fit value is the solid red line, with statistical uncertainties 
the dashed red line. The grey box is the statistical plus the systematic uncertainties. The fit values 
are shown in Table IIIH 



TABLE III: Fit values for H + system energy levels (in dimensionless units, atE). 



Irrep 


Lattice Size 


a t E 


&E,stat. 


&E,sys. 


x 2 /dof 


Q 


H+ 


20 3 x 256 
32 3 x 256 


0.291501 
0.290001 


0.000457 
0.000804 


+0.000099 
-0.000268 
+0.000418 
-0.000001 


1.003 
0.850 


0.460 
0.708 



Using the values in Table IIHI one can determine that the difference in energy between 
the two volumes is SE H + = 0.00150 ± 0.00105, where the statistical and systematic errors 
have been combined in quadrature. This sub-percent level difference is indicative of very 
small volume effects in the calculation. Also, if the most accurate data from the 20 3 x 256 
lattices is used to naively set the scale, the resulting spatial lattice spacing would be a s = 



0.12038 ± 0.00022 fm, a percent-level difference from the scale set in Ref. |28j], reflecting 
small differences from the physical point extrapolation. 
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B. s = |: Gf Irrep 

On the 20 3 x 256 lattice size calculations were also performed on the Gf irrep, which 
has only one embedding. This is an excited state of the Q~ particle with S = |, with the 
effective mass shown in Fig. and the results of fitting the plateau in Table IIVI The results 
in Fig. clearly show worse signal to noise behavior than for the H + state. This behavior 
precluded the examination of the two-G+ system. The ratio of the extracted Gf mass to 



the H + mass lattices compares within error to that extracted in Ref. [16 
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FIG. 5: (Color online) Effective mass plot for the Gf (S=^) $7 baryon calculated using 20 3 x 256 
lattice. The fit value is the solid red line, with statistical uncertainties the dashed red line. The 
grey box is the statistical plus the systematic uncertainties. The fit value is tabulated in Table ITVT . 



TABLE IV: Fit values for Gf system energy levels (in dimensionless units, atE). 



Irrep 


Lattice Size 


a t E 


&E,stat. 


&E,sys. 


X 2 /d°f 


Q 




20 3 x 256 


0.422541 


0.003754 


+0.005010 
-0.002036 


0.409 


0.931 



VI. STRANGENESS -6 SYSTEMS 

A. s = (| <g) |): The Af, E + , and T 2 + Irreps 

Two strangeness -3 H + baryons (two Q~ particles) can combine to make a strangeness 
-6 system. By forcing this system to be in a relative s-wave the angular momentum of the 
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resulting state will be entirely determined by the spin combinations allowed, which for the 
two Q~ system are the S = (the irrep) and the S = 2 (the E + and T 2 irreps) angular 
momentum states. 

For the irrep each of the H + baryons can be put into two embeddings at both 
the source and the sink, allowing for six embedding combinations: Af. u n , Af. n 22 , Af. 12 n , 

-^1-12,22) ^i^22,ii' an< ^ -^1-22,22 where Tij^i has source embeddings % and j with sink embeddings 
k and / for irrep V. Through a similar analysis of each embedding combination as was 
performed for the H + , the combinations Af. 1222 an d A\. 22 22 wer e observed to plateau well 
above the common ground state that the other combinations found. Note that these two 
are the only combinations where it is impossible to avoid a contribution from contractions 
similar to those found in the H 22 embedding combination, and thus they likely suffer from a 
similar set of excited state contaminations. Given these observations, these two embedding 
combinations are excluded from the analysis of the A^ system. The remaining embedding 
combinations are summed and result in the effective mass plots in Fig.[7]for both the 20 3 x 256 
and the 32 3 x 256 lattices, with the fit values and energy shifts (AE as defined in Eq. [3]) 
given in Table |Vj 
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FIG. 6: (Color online) Effective mass plots for the (S=0) two Q~ system calculated using (a) 
20 3 x 256 and (b) 32 3 x 256 lattices. The fit value is the solid red line, with statistical uncertainties 
the dashed red line. The grey box is the statistical plus the systematic uncertainties. The fit values 
are shown in Table [V] 



TABLE V: Fit values and Energy Shifts for Af system energy levels (in dimensionless units, atE). 



Irrep 


Lattice Size 


a t E 


&E,stat. 


PE,sys. 


X 2 /dof 


Q 


a t AE 


&AE,stat. 


At 


20 3 x 256 
32 3 x 256 


0.586235 
0.583224 


0.000843 
0.002002 


+0.000091 
-0.000348 
+0.000577 
-0.000680 


1.105 
1.086 


0.327 
0.350 


0.00323 
0.00322 


0.00124 
0.00257 



The E + and T 2 irreps correspond to the S = 2 two f2 - system. Due to the limited 
number of quark spin states available to make this spin structure, the embedding combi- 
nations that will produce a non-zero result are much more limited than in the A^ case. 
Specifically, placing either both of the source or both of the sink baryons into the first H + 
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embedding is forbidden. Furthermore, the more complicated structure inherent to this spin 
state significantly reduced the signal to noise ratio in many of the remaining embedding 
combinations, leaving usable signals only the combinations E J [ 222 and T^-vin- The effective 
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FIG. 7: (Color online) Effective mass plots of the (a) E + and (b) irreps (the S = 2 two f2~ 
system) calculated using the 20 3 x 256 lattice. The fit value is the solid red line, with statistical 
uncertainties the dashed red line. The grey box is the statistical plus the systematic uncertainties. 
The fit values are tabulated in Table EO 

mass plots for these embeddings are shown in Fig. [7] and the fit values and energy shifts are 
in Table IVII From Table I VI I one can see that the two irreps achieve statistically separate 
lowest energy states, despite coupling to states with the same set of quantum numbers. This 
would indicate that at least one, and possibly both irreps, are failing to achieve the correct 
ground state of the S = 2 two Q~ system. In both cases, however, the states achieved are at 
a significantly higher energy level than for the 5 = case, implying a much more repulsive 
channel, as expected from Pauli exclusion arguments. 



TABLE VI: Fit values and Energy Shifts for the S = 2 two Q system energy levels (in dimen- 
sionless units, atE). 



Irrep 


Lattice Size 


a t E 


&E,stat. 


&E,sys. 


X 2 /dof 


Q 


a t AE 


O 'AE, stat. 




20 3 x 256 


0.642961 


0.007136 


+0.002502 
-0.005120 


0.925 


0.514 


0.05996 


0.00719 


E+ 


20 3 x 256 


0.67256 


0.00293 


+0.00013 
-0.00329 


0.500 


0.916 


0.08956 


0.00307 



B. S = (| <g> |): The T+ and E + Irreps 

The final system examined in the calculation was the interaction between two strangeness 
-3 baryons, one in the H + irrep and one in the Gf irrep, on the 20 3 x 256 lattices. This 
system is unique in that it involves explicitly placing one baryon in in the S = \ excited 
state. The final spin combinations allowed for this combination are S — 1, which will fall 
into the irrep, and 5 = 2 which will again fall into either the E + or T 2 + irreps. The Gf 
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irrep has only one embedding, so the possible embedding combinations are determined solely 
by the remaining H + irrep baryon. In the S = 1, 7\ + case usable signals were recovered 
for embedding combinations with the source H + in the first embedding and the sink in the 
second. For the S = 2 case a signal was uncovered only for the E + irrep again with the 
source H + in the first embedding and the sink in the second. The effective mass plots for 
each are shown in Fig. [S] and the fit results and energy shifts are in Table IVHI 
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FIG. 8: (Color online) Effective mass plots of the (| <g> ^) -coupled two O system calculated using 
the 20 3 x 256 lattice for (a) (S=l) and (b) E + (S=2) irreps. The fit value is the solid red line, 
with statistical uncertainties the dashed red line. The grey box is the statistical plus the systematic 
uncertainties. The fit values are tabulated in Table EED 



TABLE VII: Fit values and Energy Shifts for the Gf <S> H + system energy levels (in dimensionless 
units, atE). 



Irrep 


Lattice Size 


a t E 


CTE,stat. 


PE,sys. 


X 2 /dof 


Q 


a t AE 


&AE,stat. 


T+ 


20 3 x 256 


0.679179 


0.002773 


+0.001087 
-0.000992 


0.389 


0.961 


-0.03486 


0.00469 


E+ 


20 3 x 256 


0.695768 


0.007049 


+0.007353 
-0.000774 


0.747 


0.650 


-0.01827 


0.00800 



VII. SCATTERING AND k cot<5 

With energy levels of each baryon and system of baryons determined in Sect. |V] and 
Sect. IVI| respectively, one can now extract scattering information from this data following 
the discussion from Sect. IHIl Returning first to the (| Cg> |) S = two Q~ system, the data 
from two different volumes will allow for two applications of Eq. [Hand, i n combination with 
Eq. El an extraction of the scattering length a. In principle the range parameter r will also be 
extracted, however that term in Eq. El will also be contaminated by contributions from all of 
the higher order terms in the expansion, and as such will be unreliable. This determination 
of the scattering characteristics of the two Q system will allow for definitive statements to 
be made on the form of interaction between these baryons in light of the conflicting claims 



of Ref. 251 and Ref. 26 . 
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Using the data from Table HTT1 and Table IVl along with Eq. [3]and Eq. [4]one can determine 
the k 2 and kcotS(k) values for the S = two Q~ system. The (dimensionful) results are 
shown in Fig. [9(a) along with the systematic and statistical errors. To obtain the scattering 
length with a correct propagation of errors, a distribution of the parameter a is generated. 
To accomplish this a series of 10,000 pairs of random values were taken from the distributions 
of k 2 for both the 20 3 x 256 and the 32 3 x 256 lattice data. These distributions were normal 
distributions defined from the mean value and statistical plus systematic uncertainty of k 2 
obtained from Table IIII| Table [V] and Eq. [3j Each random pair was then used to generate a 
pair of kcot5 values using Eq. HI which was then fit to Eq. [6] to produce a single value for the 
scattering length a. The distribution of these values is shown in Fig. Efb), with a resulting 
two Q~ scattering length in the S=0 channel of 

a g2 = 0.16 ±0.22 fm. (18) 

Note that the distribution fit to the data is that of a Lorentz distribution and not a normal 
distribution, due to the specific form that Eq. Uptakes. If one used standard error propagation 
techniques to determine the uncertainty in the scattering length, the form of the distribution 
as a Lorentz distribution rather than a normal distribution would result in the quotation of 
too large an uncertainty. 




0.0 0.5 1.0 1.5 2.0 -2 -1 1 2 



k 2 (fm" 2 ) a? = n o (fm) 

(a) (b) 

FIG. 9: (Color online) Plot of (a) kcotS and (b) the distribution of scattering lengths a for the 
S = two f2~ system. The kcot5 plot consists of the central value (black dot), statistical error 
(thick inner red line), and statistical plus systematic error (thin outer blue line) for the 20 3 x 256 
lattices (top line) and the 32 3 x 256 lattices (bottom line). The scattering length distribution 
has the probability density histogram overlaid with a Lorentz distribution fit to the data. This 
distribution has a central peak at 0.16 fm with a 68% confidence interval (equivalent to la) of 0.22 
fm. 

From Fig. [9(a) one can see that the central value of k 2 observed does not change appre- 
ciably between the two different volumes. Also, in Fig. [9(b) the distribution of the extracted 
scattering length is strongly peaked at very small values. Both of these pieces of information 
are indicative of a very weakly repulsive system. Indeed, if one operates with the assumption 
of natural sizes for the range and higher order parameters in the effective range expansion 
then the Lorentz distribution for the scattering length would provide an 79.5% chance that 
the system is repulsive and a 20.5% chance that it is attractive. Additionally, the AE values 
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in Table |V] are positive (repulsive) and small within the la error band for both 20 3 x 256 
and 32 3 x 256 lattices. Thus, from our current lattice calculations, we find evidence that 
the system is consistent with the weakly repulsive scenario in Ref. M 

and inconsistent with 

the deeply bound state found in Ref. [25[. Ultimately, more calculations are required to 
acquire a full error budget of the systematic, but these systematics for the QQ system are 
not expected to be appreciable for the reasons mentioned earlier. 

Finally, for the other two baryon systems studied on the 20 3 x 256 lattices a calculation 
of k 2 and kcotS can be made in a manner similar to that for the 5 = two Q~ system 
detailed above. However, as these calculations were performed on only one lattice size, a 
reliable scattering length cannot be obtained and one can only quote the inverse of kcotS as 
a proxy for the scattering length. The inverse of kcotS is plotted for each of the two baryon 
systems studied on the 20 3 x 256 lattices in Fig. [10J The S = two Q~ system already 
discussed is the middle of the five points shown. 
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FIG. 10: (Color online) Combined results for kcot5 for all two baryon systems studied on the 
20 3 x 256 lattices. These systems are (from left to right) the S = 1 Gf H + system (orange circle), 
the S = 2 G J [H + system (red square), the 5 = H + H + system (brown diamond), the S = 2 
H + H + system in the irrep (pink downward triangle), and the S = 2 H + H + system in the E + 
irrep (purple upward triangle). The errors shown are statistical only. 



Of note in Fig. [TO] is that the GiH + excited state system appears to be attractive in 
both the S = 1 and the S = 2 channels. However, without additional lattice volumes one 
cannot make a claim as to whether these states would be bound or simply scatter attractively 
(resonances). Also, in nature the state would decay strongly to a H baryon and a kaon 
in a relative p-wave with the release of several hundred MeV of energy (assuming a Gf mass 
of 2500 MeV[l6[), a far lower energy arrangement than any putative GfH + bound state. 
Also, the ordering of the GfH + systems does seem to follow the general trend observed and 
expected in the other two Q~ system calculated where the lower spin state results in a more 
attractive/less repulsive interaction. 

For the S = 2 H + H + system in Fig. [TU1 the two different irrep calculations clearly 
obtain two different states. It is notable that the two states lie on different branches of 
the kco\8 function, and that neither are on the principle branch that the S = states 
are on. This would indicate both that the S = 2 states are strongly repulsive (which one 
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could qualitatively predict from the Pauli exclusion principle) and that it is possible that 
neither of the extracted states are actually the S — 2 ground state. If the interaction is 
repulsive enough then the volume in which the objects are contained may not be able to 
physically contain the two baryon system, leading the scattering interaction to be pushed to 
higher branches of kcotS. In a regime where the interaction fits within the lattice volume, 
the interaction energy should scale inversely with the volume, with any stronger scaling 
indicating a volume that is too small. A calculation of the S = 2 H + H + system at a second 
volume would point toward a resolution of this question. 



VIII. CONCLUSION 

We have performed the initial s-wave lattice calculations of the two spin-3/2 Q~ system 
for S—0 and 2 channels, as well as spin- ^ /spin- 1 coupled two-f2 system for S=l and 2 
channels. Using lattice configurations with pion mass m n ~ 390 MeV, calculations were 
performed in a volume with L=2.5 fm, and in the S = case, also with L=3.9 fm. The S=0 
results demonstrate that the two-fi system in this channel is most likely weakly repulsive, 
with a scattering length 

= 0.16 ± 0.22 fm. (19) 

Also, the energy of interaction was positive and small within errors for both L=2.5 fm and 
L=3.9 fm, whose values are states in Table IV1 As our calculations rely on no phenomenology, 
we assert that these findings should provide significant evidence supporting the conclusion of 
a weakly repulsive system in Ref. j26j, as opposed to a deeply bound state in Ref. jHf. These 
results also provide an interesting complement to previous studies EMI, |30j of hyperon 



interactions, where many of the interactions have been found to be attractive and contain 
bound states at a pion mass of 390 MeV. The difference between the evidence for other 
bound hyperon states and the conclusion in this work of a weakly repulsive QQ state may 
simply reflect a much stronger influence of light-quark dynamics in the valence sector of the 
AA and systems. Further studies at different pion masses approaching the physical point 
are needed to gain a better understanding of the similarities and differences of each of these 
systems. 

Central to this work was the use of interpolating operators of definite lattice cubic sym- 
metries, which allowed us to look at exotic channels of the two-f2 system. In some of these 
channels we find evidence for an attractive interaction. Future studies with larger statistics 
and multiple volumes will better elucidate the behavior of two-f2 systems in these channels. 
The minimal light-quark dependence of the Q system implies that the results presented in 
this work should be 'near-physical'. Extension of this work to configurations with lighter 
pion mass, as well as the development of an xEFT for this system, will quantify this state- 
ment. 
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